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Abstract: Piezoelectric material has been emerging as a popular building block in MEMS 
devices owing to its unique mechanical and electrical material properties. However, the 
reliability of MEMS devices under buckling deformation environments remains elusive 
and needs to be further explored. Based on the Talreja's tensor valued internal state 
damage variables as well as the Helmhotlz free energy of piezoelectric material, a 
constitutive model of piezoelectric materials with damage is presented. The Kachanvo 
damage evolution law under in-plane compressive loads is employed. The model is applied 
to the specific case of the postbuckling analysis of the piezoelectric plate with damage. 
Then, adopting von Karman's plate theory, the nonlinear governing equations of the 
piezoelectric plates with initial geometric deflection including damage effects under 
in-plane compressive loads are established. By using the finite difference method and the 
Newmark scheme, the damage evolution for damage accumulation is developed and the 
finite difference procedure for postbuckling equilibrium path is simultaneously employed. 
Numerical results show the postbuckling behaviors of initial flat and deflected 
piezoelectric plates with damage or no damage under different sets of electrical loading 
conditions. The effects of applied voltage, aspect ratio of plate, thick-span ratio of plate, 
damage as well as initial geometric deflections on the postbuckling behaviors of the 
piezoelectric plate are discussed. 
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1. Introduction 

The use of piezoelectric materials in intelligent structures has received considerable attention in 
recent years due to the intrinsic direct and converse piezoelectric effects. Piezoelectric materials have 
been used as sensors or actuators for the control of the active shape or vibration of structures. Defects 
such as microcracks, voids, dislocations and delamination are introduced in piezoelectric materials 
during the manufacturing and poling process. The existence of these defects greatly affects 
the electric, dielectric, elastic, mechanical and piezoelectric properties of the piezoelectric materials, 
especially the service life of piezoelectric structures. When subjected to mechanical and electrical 
loads, these defects may grow in size and cracks may propagate leading to premature mechanical or 
electrical fatigue failure. Therefore, it is important to understand the growth of these defects, the 
damage accumulation and the overall effect of these defects on the average mechanical and electrical 
properties of piezoelectric structures. 

Damage in fiber-reinforced composite materials has been extensively investigated, and many 
theories have been established and used to predict the life of composite structures. Based on the 
framework of irreversible thermodynamics with internal state variables, Talreja [1] developed a 
phenomenological theory for composite laminated plates. In his study, the Helmholtz free energy was 
expanded into a polynomial in terms of elastic strains and damage variables to obtain the stiffness-damage 
relations. Utilizing a continuum mechanics approach, Allen et al. [2,3] developed a model for 
predicting the thermomechanical constitution of initially elastic composites subjected to both 
monotonic and cyclic fatigue loading. Valliappan et al. [4] established the elastic constitutive 
equations for anisotropic damage mechanics, and the implementation of these constitutive equations in 
the finite element analysis was explained. By defining damage variables as the material stiffness 
reduction, Ladeveze and Dantec [5] formulated the constitutive equations and the corresponding 
damage evolution laws of the elementary ply for laminated composites that can be used to describe the 
matrix micro-cracking and fiber/matrix debonding. Schapery and Sicking [6] discussed the 
homogenized constitutive equations for the mechanical behavior of unidirectional fiber composites 
with growing damage, and the emphasis was on resin matrices reinforced with high modulus elastic 
fiber. Zhang et al. [7] investigated a computational model for the damage evolution of engineering 
materials under dynamic loading, and two models for dynamic damage evolution of materials in 
general anisotropic damage state were presented. Moore and Dillard [8] have observed time dependent 
growth of transverse cracks in graphite/epoxy and Kevlar/epoxy cross-ply laminated at room 
temperature. Luo and Daniel [9] have shown that the macroscopic mechanical behavior of 
unidirectional fiber-reinforced brittle matrix composites can be correlated explicitly with the 
microscopic deformation and damage. 

Modeling and analysis of multilayer piezoelectric beams and plates have reached a relative maturity 
as attested by the numerous papers. Mindlin [10] presented the theory of piezoelectric crystals plate 
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considering shear and bending. Chandrashekhara and Tiersten [11] developed general piezoelectric 
nonlinear theory and detailed the vibration equations of different piezoelectric crystals. 
Chandrashkhara, Tenneti [12] and Zhou et al. [13] investigated the dynamic control of laminated 
piezoelectric plates by the finite element (FE) method. Wang and Rogers [14] presented a model for 
laminated plates with spatially distributed piezoelectric patches. Tzou and Gadre [15] analyzed thin 
laminates coupled with shell actuators for distributed vibration control. Xu et al. [16] analyzed the free 
vibration of laminated piezothermoelectric plate based on the 3D theory. Mitchell and Reddy [17] 
proposed the theory of the laminated piezoelectric plates by using classical plate theory and simple 
third-order theory, respectively. Noor and Peters [18] presented postbuckling analysis of multilayered 
composite plates subjected to combined axial and thermal loads and investigated the effects of 
mechanical and thermal loading on postbuckling behaviors of composite plates. Based on transverse 
shear-deformable theory, Librescu and Souza [19] gave the postbuckling analysis of geometrically 
imperfect flat panels under combined thermal and compressive edge loadings. Shen [20] investigated 
the postbuckling behaviors of laminated plates with piezoelectric actuators under complex loading 
conditions based on Reddy 's higher order shear deformation plate theory. Oh et al. [21] studied 
thermal postbuckling behavior of laminated plates with top and/or bottom actuators subjected to 
thermal and electrical loads. In their analysis all the static/dynamic behaviors of laminated plates were 
investigated without considering the damage effects which would greatly influence the mechanical 
behaviors of smart structures. To the best of the author's knowledge, up to now postbuckling analysis 
of piezoelectric structures considering the damage effects has rarely been investigated and reported. 
For example, Aydin [22] studied the dynamic characteristics of functionally graded beams with open 
edge cracks, in which an analytical method was proposed to determine the free vibration of beams with 
any number of cracks. Mao et al. [23] studied the creep buckling and post-buckling of laminated 
piezoelectric viscoelastic functionally graded material (FGM) plates by adopting the Boltzmann 
superposition principle, and the nonlinear creep buckling governing equations of the laminated 
piezoelectric viscoelastic FGM plates with initial deflection were derived on the basis of the Reddy's 
higher-order shear deformation plate theory. Hamed [24] presented a nonlinear theoretical model for 
their bending and creep buckling analysis. The model accounted for the viscoelasticity of the materials 
using differential-type constitutive relations that were based on the linear Boltzmann' s principle of 
superposition. Xu et al. [25] presented a high performance and simple structure bi-stable piezoelectric 
energy harvester based on simply supported piezoelectric buckled beam. Cottone et al. [26] 
investigated an approach for piezoelectric beams by exerting an increasing axial compression and 
demonstrated that the numerical model and experimental results were in good qualitative agreement. A 
constitutive model of fully coupled electro-magneto-thermo-elastic multiphase composites has been 
proposed by Aboudi [27]. In his works, the linear displacement, electric potential and magnetic 
potential are adopted, which can't predict the micro fields precisely. Bansal and Pindera [28] proposed 
a unified macro-and micro -mechanics failure model with method of cells as finite-volume direct 
averaging micromechanics (FVDAM). Based on the FVDAM theory, Sun et al. [29] built a unified 
macro- and micro-mechanics constitutive model of fully coupled fields in composite materials by 
high-order displacement, electric potential and magnetic potential. 

In the present study, a new constitutive model for piezoelectric materials using the Talreja's tensor 
valued internal state damage variables and the Helmhotlz free energy of piezoelectric material is 
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presented. This model is then applied to a specific case of postbuckling analysis of piezoelectric plates 
under in-plane compressive loads. By adopting von Karman's plate theory and using the finite 
difference and the Newmark scheme, the damage evolution for damage accumulation is developed and 
the finite difference procedure for postbuckling equilibrium path is simultaneously employed. In the 
numerical examples, the effects of variation in the load parameters, damage influences and geometric 
parameters of the plate on postbuckling equilibrium paths are discussed. 

2. Basic Equations 

2.1. Constitutive Equations for Damaged Piezoelectric Materials 

Consider a representative volume element of a piezoelectric solid with a multitude of damage 
entities in the form of microcracks, as shown in Figure 1. As discussed in Talreja, two vectors are 
needed to define each damage entity. These are the damage influence vector a, and the n t unit normal 

to the damage entity surface. The damage influence vector represents an appropriately chosen effect of 
the damage entity on the surrounding medium. With these two vectors, a damage entity tensor d tj is 

formed by taking an integral of the diad a i n . over the surface of the damaged entity: 

d ij=\ s ^j dS (1) 

where S is the damage entity surface. 
Figure 1. Representative volume element with internal damage variables for piezoelectric materials. 




Now if there are n distinct damage modes in the representative volume element (e.g., intralaminar 
cracks in different orientations, etc.) denoted by k =1,2,... n, a damage tensor can be defined for each 
mode as: 

^4z^ (2) 

where V r is the volume of the representative volume element and 9 k represents the number of damage 
entities in the kth damage mode. The tensor co.j is an unsymmetrical tensor in general. However, we 
can represent the vector a i along the normal and tangential directions at any point on the surface of the 
damage entity and write: 

d ij= d l+ d l (3) 
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where d\ } =^an i n } dS and d 2 . = £ bm i rijdS .in which, a and b are the magnitudes of the normal and 

tangential projections of vector a i respectively, and vectors n i and m i are unit normal and tangential 
vectors, respectively. Thus the damage tensor^, can be written as: 



at=af+atf (4 ) 



where cof =^i\ /V r and cof =X( d ?A /V r- 

Physically, the damage tensor co lk represents the effects of crack opening on the surrounding medium 
whereas the damage tensor co 2k represents the effects of sliding between the two crack faces. In many 
situations, the sliding between the crack faces can be negligible, e.g., for intralaminar cracks constrained 
by stiff plies, and hence we assume co 2k = 0 . This implies co k = co) k which is a symmetric tensor. 

For the case of damaged piezoelectric material without temperature effect where the damage is 
represented by internal state variables, the Helmholtz free energy of piezoelectric material can be 
written as a function of the transformed elastic strains, the electric field vector and damage internal 
variables, that is: 

H=H(e il ,E i ,a>*) (5) 

The transformed stress components cr and the electric displacement components G i at any fixed 
damage state are now given by: 

a = — 

d £ij 

8H(e u ,E„^) (6) 

When the damage induced by the cracks in the piezoelectric material has the ortho tropic property, 
the irreducible integrity bases for a scalar polynomial function of two symmetric second rank tensors 
can be expressed as [30]: 

2 2 2 

1 ' ^22 ' ^33 ' ^23 '^31 '^12 ' ^1 2 ^23 ^3 1 ' 

^^k ^ k ^^k f^^k\2 /~~k\2 / ~~k\2 ~~k^^k^^k 

co n , co 22 , co 33 , {0) 23 ) , {co 3l ) , yCO\ 2 ) > (0 Yi (0 n (0 z\ 

£ 23 C0 23 , S 3l C0 3i , S 12 CO l2 , C0 23 S l2 S l3 , C0 3l S 32 S l2 , CO l2 S l3 S 23 , ( 1'^' ^) (7) 

k k k k k k 

£ 23 CO l2 CO l3 , S 3l CQ 32 CQ l2 , S l2 CO l3 C0 23 

E\ , E 2 , E 3 

where n is the number of the cracks' direction in the material. For a piezoelectric single-layer plate, 
the local coordinate system o — 1 — 2 — 3 is selected, in which 1, 2 denote the two principal direction of 
the piezoelectric plate, 3 is vertical to the midsurface. According to the Kirchhoff hypothesis for plate 
s u =£ 23 =0 and applying Voigt notation to describe strains and damage variables, the bases of 
invariants can be further written as: 

s v e 2 , e 3 , s 2 , co k , co 2 k , co k , (co k ) 2 , (co k ) 2 , (co k f, 

co k co k co k , £ 6 co k , £ 6 co k co k , (k = 1,2, n) (8) 

E-y, E 2 , E 3 
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Using the above stated irreducible integrity bases, the Helmholtz free energy of piezoelectric 
materials can be expressed as a quadratic expression of the strains or the electric field intensity, a 
mixture quadratic expression of strains and electric field intensity and a linear expression of damage 
variables [31] as follows: 

H — C\£\ + C 2 g 1 g 2 + C 3 £ 2 + C^Sfr + C?s^ + C^s t s 3 + CjS 2 s 3 ~ + ^2^2 ~^ 
k 3 E 3 + k\E x E 2 + k°E 2 E 3 + *r 6 °£ 3 £i) - (e\E x s x + e\E x e 2 + elE x s 3 + e°E 2 s l + 



6 ^ Ei 2 £ '2 £ ^ ^^2 £ 3 ^7 ^'3 £ \ ^8 ^^3^*2 ^"9 ^^3 £ 3 ^ ^^"l ^* 1 £ \ ^^2 ^1 ^^3 

<fc=l 

. s~ik J2 k , s~ik Jl k . /-ik Jl ,,ifc . /-ik Jl , s~\k Jl k . /~\k Jl . /-ik Jl ^k , 

-\-C 4 S 2 CO l +t 5 £ 2 C0 2 +t 6 £ 2 C0 3 -\-C 1 S 3 CO l +t H £ 3 C0 2 + t 9 £ 3 C0 3 +C |0 ^ 6 ^ 1 + 

Ck Jl ^Jt . /-ik Jl ^k , /-ik _ „ , , /-ik _ _ k , /-ik „ „ . ✓-tA' „ „ k , /-ik _ „ , £ 

11^6^2 + ^12^6^3 +^i3^ 1 ^ 2 <y i + ^14^2^2 + ^15^2^3 + ^16^2^3^ + ^\1 £ 2 £ Z G) 1 

+ ^18 £ 2 £ 3^3 + ^19^3^1 "'l + ^O^^l"^ + ^21^3^1^3 + ^22^6^1^6 + ^ y 23 £ b £ 2 ( ^6 

C^^o* ) - ^ (k\E\(q\ + k 2 E?co 2 + k* El col + K^l 0 ^ + *5 ^2 ^2 + K 6 E l W l + 

k=\ 

kIeIco\ + kIe]co\ + k\eIcoI + k\^E 2 co\ + k\ x E x E 2 co\ + k\ 2 E x E 2 coI + k1 3 E 2 E 3 &1 

n 

+ Ki 4 E 2 E 3 co k 2 + k1 5 E 2 E 3 ojI + K^E^col + k^E^E^I + k\ % E 3 E x (oI ) - ^ {e\E x sjo\ 



(9) 



+ <?2-E'i£'i<^2 +^3£ , i^i<«3 +e4£ , 1 £ 2 <y* +elE 1 e 2 a>l +e l ^E l £ 2 6)l +e k 1 E x e 3 o)\ +elE l £ 3 o; 



~\~ @q Ei^£ 3 C0 3 ~\~ (?^q£j 2 £^C0^ €-^^Ei/2£^(Q 2 ~\~ £y2Ei 2 £^X) 3 H - £y 3 Ei 2 £ 2 Q)^ ~\~ (?^^£j 2 £ 2 C0 2 ~\~ &y^Ei 2 £ 2 C0 3 



~\~ &y^Ei 2 £ 3 Q)^ @^Et 2 £ 3 C0 2 ~\~ Gy^Ei 2 £ 3 Ct) 3 ~\~ £ 3 £ ^CO^ ^2,0^^3 £ \^^l ^i\^ i 3 £ \^^3 ^2.1^'3 £ 1^^\ 

+ e 23 E 3 £ 2 a> 2 + e 24 E 3 £ 2 a> 3 +e 25 E 3 s 3 0J 1 +e 26 E 3 s 3 Q) 2 + e 21 E 3 e 3 co 3 + e 2S E l s 6 a> 6 +e 29 E 2 £ 6 o t 



k 

6^6 



where C°(i = 1,2, -"7) are the elastic material constants without damage, = l,2,---9) are the 
piezoelectric constants without damage, icf(i =1,2, •••6) are the permittivity matrix constants without 
damage, C.(i = 1,2, •••24) are the material constants with damage; ef (/ = l,2,---30) are the 
piezoelectric constants with damage, k\ (i = 1,2,---18) are the permittivity matrix constants with 
damage,/? is the density of piezoelectric material, P 0 is a constant, Pj is a linear function of strains, P 2 is a 
linear function of damage variables andP 3 is a linear function of the electric field intensity. Then the 
stresses and the electric displacements can be expressed as: 

dH 



ap de. 



= [C° +Yc k ]e -[e° +Ye k ]E 

L pq / j pq 1 a L pm / \ pm s m 



G_=- 



p 

dH 



k=l 



k=l 



(10) 



k=l 



k=l 



where C° pa , C k pg , tc^ and/c^ are all symmetric matrixes having the forms as follows: 

"2c, 0 c° c° 0 



[C° ] : 

pq 



2C° 3 



2C 5 ° 



0 
0 

2c: 



(11) 
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2C,X+2CX + 2C 3 X 
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(12) 
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(13) 



(14) 



e\ co\ + e k 2 co 2 + e\ co\ 
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1 

e 



co k +e 23 co 2 +(\,oj. 



-24^3 



e k 0 co k 



k k 



k k k k 

+ K n C0 2 + K n G>l 



in '-11— 2 " 

2k\ co\ + 2k\(q\ + 2k\(o\ 



h k k k 

+ K\ 1 C0 2 +K" 18 ^3 
K k n C0\ + K k 4 C0 2 + K k 5 CO k 

2K k co\ +2K k co k +2K k co k 



(15) 



(16) 



Assuming that there is only one damage mode in the representative volume element, the relations of 
the strains, the stresses, the electric field intensity and the electric displacements in Equation (10) can 
be simplified as: 



a = [C° +C 1 ]e -[e° +e ]E 

P P<Z PQ <? P m P m m 

m [^qm ^qm -I ^ q *-^mn ^mn 



(17) 



where Cj 9 ,eJ m and< n are the same as before. C k pq , e pm , K k mn are replaced by C x pq 

In present study, consider that the matrix cracks in the piezoelectric plate are parallel to the 



Equations (12), (15) and (16) can be simplified as: 



2Cj »[ C l3 a> l C\ 9 CQ\ 
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2k 1 co 1 



(18) 



(19) 



(20) 
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Due to the fact the cracks are parallel to the coordinate plane 2 - 3 , the effect of the damage on the 
stiffness of the plate in this coordinate plane2-3can be neglected, which means the component C l n 
and C^of stiffness matrix due to damage effect are negligible. Then matrix (18) can be further 
simplified as: 



2CjO, CjjOj C 19 Oj 0 
0 0 
0 0 
2C 10 « 1 



(21) 



Letting cr 3 =0 based on plane-stress assumption and using Equation (17), the constitutive relations 
with damage of the piezoelectric plate for the plane stress problems are obtained as follows: 
[C ]s -[e ]E =[C° +C 1 ]s -[e° + e l ]E 

L na J a L nm J m L no no J a L nm nm J m 



cr 



D =[e ] T e +[k ]E =[e° +e l fe +[k q +K 1 ]E (p,q = 1,2,6 m,n = \,2,3) 

m L qm- 1 q L mn J n L qm qm J q L mn mn J n \t"> x l ' ' ' ' ' / 



(22) 



where: 



id 
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(28) 



In the present research, the Kachanvo damage evolution law [32] is adopted for an arbitrary point i 
of the piezoelectric plate with damage: 



do){ 
dt 



l -co! 1 . 1 

W q \<ju(l-col)(j f 



(29) 



0 



where B , m and ju are the material constants, <j eq is an equivalent stress which is based on certain failure 
criterion, a f is the limit stress. 

The relations between the electric fields E X ,E , E z and the electric potential cj> in the Cartesian 
coordinate system are defined by: 

E x=-ts> E y=-0,y> E z= -</>,;. (30) 

For the piezoelectric plate, only thickness direction electric field E v is dominant. If the voltage 
applied to the piezoelectric plate with piezoelectric effect in the thickness only, then: 

E z = v J h (31) 
where V is the applied voltage across the thickness of piezoelectric plates. 



2.2. Basic Equations of Piezoelectric Plates 

Now, consider a thin piezoelectric plate with transverse cracks subjected to uniformly distributed 
in-plane compressive load P in the ^-direction, having thickness h , length a in the jc-direction, width 
b in the v-direction shown in Figure 2. The reference surface defined by z = 0 is set on the middle 
surface of the undeformed plate. 
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Figure 2. Geometric configuration of a piezoelectric plate with transverse cracks under the 
uniform compressive in-plane loads. 
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cracks 
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cracks 




Setting u ,v and was the displacement components of an arbitrary point on the mid-surface along the 
direction of x , y and z , respectively, and denote w as the initial geometric deflection. According to 
classical nonlinear theory, the strain components s° , s^mdy^ of the mid-surface [33] can be written as: 



S l = U , x + \ W , 2 + W ,x W ,x> £ y = V ,y + \ W ,y + W ,y W ,y> rly = U y + V , + W , W y + W ,x W ,y + W ,y W ,x 

and the curvatures k x , x y and K xy of the mid-surface as: 

K x = ~ W ,xx . K y = - W ,yy • ^xy = 

then the nonlinear strain-displacement relations are expressed as follows: 

e x =e° x +ZK x , £ y =s° y +zrc y , y xy =rl+zrc xy 



(32) 



(33) 



(34) 



Suppose the damage variable remains constant through the thickness of plate. Denoting N x ,N y , N xy as 
the membrane stress resultants and M x ,M y ,M xy as the stress couples of the plate, according to the 
classical nonlinear plate theory, the nonlinear governing equations of the piezoelectric plate with initial 
geometric deflection [33] can be written as: 



N +N =0 

x,x xy,y 

N +N =0 

xy,x y,y 



(35) 



M +2M +M +N O +w ) + 2N (w +w ) + N(w +w )=0 

x,xx xy,xy y,yv x v ,xx ,xx / xy ^ ,xy ,xy / y v ,yy ,yy 7 

Using Equations (22) and (34), the following constitutive equations can be obtained: 
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(36) 



(37) 



where N X ,N X ,N X . and M x ,M p ,M p represent the component resultants and couples due to the 
piezoelectric effect, respectively. The stiffness coefficients A tj and Z> of the piezoelectric plate are 
defined as follows: 
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\=&A> D iJ=\!y/ C iA (U = 12,6) (38 ) 
The resultants and couples due to the piezoelectric effect can be written as: 



ATf = f /2 e 3l — dz , N% = \ e 32 ^ dz , AT = 0 

m;=o, M''=0, m?=o 



(39) 



Introducing the following dimensionless parameters: 

Z = ±r, = lu = ± V = V -, W = ™,W = ™, 
a b a b h h 

A - Au A - A>2 A - Al2 A - ^ 6 ;--!=- n - ° n 
C n /i C H /i C n /i C h /j a o C n /i 

D - D - D — C — C 

j ) _ _^12_ 7-) _ _^22_ 7-j _ _t_66_ _ *-12 _ ^11 = e» _ _ ~ / 

~~ C*h 3 ' ~~ C*A 3 ' ~~ C*A 3 ' C* ' ~~ C* ' 9 ^* ' ° s (^O) 

aw = — 4-, «i2 = — 4-, #22 = -^=-. ass = — — . e., = — . e,„ = — 
Qi Qi 




_ * ' 


(7 f = 


c" 


* > 

e 3l 


032 = 


^31 



By using Equations (22), (32-39) and (40), the dimensionless nonlinear governing equations of 
piezoelectric plate with initial geometric deflection under compressive loads in-plane including 
damage effects are obtained and expressed in terms of u,V andiv as follow: 

±V 4 +\X 2 W 4 W„ +A l A 2 W^ ri + X l X 2 W 4 W^ + A m {^U m +V 4l] +^ 2 W,W ? , +^W^W m + 

+HX 2+ W,) + 4(^„ H 2 ^,w„ +^ 2 ^„)+iA-, 6 ,,(^^ + 

+ + +f ^ + + 4 2 ^% + ( 41 ) 

I'^W, + 4 2 ^„ + A, 2 W ^ + ^ 2 W,W g ) - = 0 

*2 4 0 22m W m + 2D 22ji W^ +D 22 W m J + A 1 2 [A 11 (U 4 + ^ 2 W/ + + + W 

- ^(V^ + ^) + 2^A 66 (^,„ + £^ + + + /.x.U- U'.KU- . 

+ W^) + X 2 2 [A l2 {U 4 +\\ 2 W 4 2 +KW 4 W 4 ) + A 22 (V, 1+ \X 2 W; + A 2 2 W„WJ-^V.W Jn +W^) = 0 
Suppose the boundary of the piezoelectric plate is simply movably supported, the dimensionless 
boundary conditions can be expressed as: 

£ = 0,1: V = 0,A u {U^\\ 2 W^A x 2 W^)-7, l V e =-P (42) 
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W = 0, =0 

7 = 0,1: u = 0, A 22 (^+H 2 ^, 2 +^ 2 ^,)-4K=0 
W = 0, W m =Q 

The dimensionless damage evolution equation of the piezoelectric plate subjected to the uniformly 
in-plane compressive load can be written respectively as follows: 



dco\ 
~d7 



I -col 



0 



W eq \<ti<\-a{)5 f 



(43) 



Taking the mid-surface normal stress of the piezoelectric plate as the equivalent stress o eq that is 
parallel to the fibrous direction, it can be presented as: 



& eq =Cn (U 4 + \ \ 2 W/ + ^ Z WJV^ ) + Cn(V t] +\ A 2 Z W/ + AfWJY^ ) - e 3 \V e 



1 3 2 T 



(44) 



3. Solution Methodology 

Suppose the dimensionless initial geometric deflection is taken as: 

W = W 0 sin pnB, sin qrcr] (45) 

where p and q are the mode number in the | -direction and rj -direction of the piezoelectric plate, 
respectively. Since the load and the structure are symmetric, only one quarter of the plate needs be 
considered. So the domain of the problem is selected as0<<^<l/2,0</7<l/2. 

To seek the approximate solutions of the governing Equations (41) which satisfied the boundary 
conditions (42), the unknown functions U,V and Ware separated both for space and for time. The 
finite difference method is used for space, and the partial derivatives with respect to the space 
coordinate variables are replaced by difference form. The timer is equally divided into small time 
segments Ar, and the whole equations are iterated to seek solutions. At each step of the iteration, the 
nonlinear items in the equations and the boundary conditions are linearized. For example, at the step j , 
the nonlinear items may be transformed to: 

(x-y) J =(x) r (y) Ji ( 46) 

where (y) s is the average value of those obtained in the preceding two iterations. For the initial step of 
the iteration, it can be determined by using the quadratic extrapolation, i.e.: 

(y) Jp = A( 30,_, + B( y) ,_ 2 + C( y),_ 3 (47) 

and for the different step of the iteration, the coefficients A, B and c can be expressed as follows: 

J = \: A = 1,5 = 0,C = 0 
7 = 2: A = 2,fi = -l,C = 0 (48) 
J>3: A = 3,B = -3,C = 1 



Moreover, using the Newmark scheme, the inertia in Equation (41) can be expressed as follows: 
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At 



(W r ) y =(W r ), / _ 1 +i[(W rr ).,_ 1 +(WJ y ]Ar 

W/ = W,_, + (W r ) y _ 1 Ar + {[(W rr ),,_ 1 + (WJ,,](Ar) 2 



(49) 



For every time step, the iteration lasts until the difference of the present value and the former is 
smaller than 0.1%, then continue the calculation of the next step. 

4. Numerical Results 

4.1. Comparison Study 

To ensure the accuracy and effectiveness of the present method, a test example was calculated for 
postbuckling analysis of isotropic rectangular plate with initial geometric deflection. Comparison of 
postbuckling response curves for isotropic rectangular plate with initial geometric deflection is shown 
in Figure 3. The boundaries of the plate are clamped movable edges. W 0 denotes the center deflection 
of the plate. The close agreements between the present results and those of reference [33] demonstrate 
the present method is accurate and effective. 

Figure 3. Comparison of postbuckling response curves for isotropic rectangular plate with 
initial deflection (W 0 = 0.1, v = 1/3) . 




4.2. Parametric Study 

To study the piezo-effects and damage effects on the postbuckling behavior of the plates, several 
numerical examples were solved for initial flat and deflected plate. A piezoelectric plate consisting of 
the PZT-5A including initial damage is considered for postbuckling analysis. The material properties 
of PZT-5A are given as follows: 
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E n = E 22 = 61.0GPa, E 33 = 53.2GPa, v n = 0.35, u l3 = o 23 = 0.38 
G 12 = 22.6GPa, G l3 = G 23 = 21.1GPa, p = 7.75xl0 3 (kg/m 3 ) 

« 31 = e 32 = 7.209C/m 2 , e 33 = 15.118C/m 2 , e 15 = e 24 = 12.322C/m 2 ( 50 ) 
k u =ic 22 =1.53xl0" 8 F/m, ^ =1.5xl0" 8 F/m 

When the effect of damage is omitted and the linear strain-displacement relations are adopted, the 
dimensionless governing equation corresponding to Equation (41) is presented as: 

-\ 2 {P + e^Ve )(W% +W#) = 0 (5 1 } 

The corresponding dimensionless boundary conditions of the simply movable supported plate can 
be written as: 

| = 0, 1: W = 0, Wg=0 

77 = 0, 1: W = 0, W TO =0 (52) 

Considering a harmonic displacement solution for this buckling problem (51), the displacement that 
satisfies with the boundary conditions (52) can be expressed as: 

Wig, 77) = W Q sin mn£, sin nnr] (53) 

Substituting Equation (53) into the Equation (51), the buckling loadP m „of a perfect plate and the 
relation of the center deflection and the compressive load in-plane of the plate with initial geometric 
imperfection can be obtained, respectively, as: 

P, m =7r\\WDn+2A 2 2 n 2 Dn+^D 2 2+41 2 2 n 2 D6 6 )-e 3 ~ 1 Ve (54) 

Wii(P) jjL^m (55) 

mn 

The least critical buckling load P cr of the piezoelectric plate is determined by applying Equation (54) 
for the buckling mode (m, n) = (1, 1) , as: 

P cr =7r 2 (A l 2 Dn + 2A 2 2 Dn+^D22+4A 2 2 D6 6 )-e^ l Ve (56) 

When the geometric parameters are given as \ =0.1, Wo =0.1, the critical buckling loadP cr is 
obtained as 0.03290 and 0.03272 by using the Equation (56) and the algorithm in the present paper, 
respectively. Figure 4 shows the relations of the center deflection W Q of the piezoelectric plate and the 
in-plane compressive load P without considering the damage effect. It can be seen that the nonlinearity 
of the plate has great influence on the postbuckling paths of the piezoelectric plate with initial deflection. 
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Figure 4. Response curves of the centre deflection vs. compressive loads in plane without damage. 
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A parametric study has been carried out and typical results are shown in Figures 5-14. It should be 
appreciated that in all figures W 0 , W 0 , V e and P cr denote the dimensionless maximum deflection, the 
dimensionless maximum initial deflection of the plate, the dimensionless applied voltage acted upon 
the plate and the dimensionless least critical buckling load. In all examples the least critical buckling 
load is taken as P„ = 0.03272 , which was calculated without considering piezo-effects and with the 
geometric parameters \ = X, = 0. 1 . 

When the damage effect is in consideration, the material parameters related to damage in all 
examples are taken as: 

5=0.5 , n = \.\ , a f =4xl0" 2 , a^ = -0.1 , «^ = -0.04 , «^ = -0.1 , a^ = -0.04 , /^ = -0.1 , =-0.1 , 
jU = 0.6. 

Figure 5 shows the postbuckling response curves for an initially flat and deflected piezoelectric 
plate without damage under different electrical loads and Figure 6 shows the effects of electrical loads 
on the postbuckling response curve of a piezoelectric plate without damage under two initial deflection 
conditions, respectively. The geometric parameters are given as \ =X 2 =0.1 .Three electrical load 
conditions, referred as 1, 2 and 3, are considered. It can be seen that the negative control voltage 
results in the increase of the buckling load and the decrease of postbuckled deflection at the same 
compressive loads. In contrast, the positive control voltage decreases the buckling load and induces 
larger postbuckled deflections. It can be concluded that the positive control voltage acting upon the 
piezoelectric plate is equivalent to a compressive piezoelectric force acting in the in-plane direction of 
the plate to some certain extent, which leads to the smaller buckling loads. 
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Figure 5. Comparisons of postbuckling response curves for initially flat and deflected 
piezoelectric plate without damage under different electrical loads. 
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Figure 6. Effect of electrical loads on postbuckling response curves of piezoelectric plate 
without damage under two different initial deflections. 
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Figure 7 shows the effect of initial deflections on the postbuckling response curve of piezoelectric 
plate without damage. The geometric parameters are given as \ =A 7 =0.1. It can be seen that the larger 
the initial deflections of the plate, the larger the postbuckled deflection of the plate under the same 
compressive load, and that the postbuckled deflections of the plate under different initial deflections 
will reach the same value with the increase of postbuckling loads. 
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Figure 7. Effect of initial geometric deflections on the postbuckling response curves of 
piezoelectric plate without damage. 
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Figure 8 shows the effect of thick-span ratio of the plate on postbuckling response curves of 
piezoelectric plate without damage under different electrical loads. It can be found that the buckling 
loads increase with the increase of the thick-span ratio of plate, and that the control voltage has a small 
effect on the postbuckling behaviors of the plate with lower thick-span ratio. 

Figure 8. Effect of thick-span ratio on postbuckling response curves of piezoelectric plate 
without damage under different electrical loads (the inset figure is a zoom-in snapshot of 
the region around orgin point to depict the difference of three cases). 
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Figure 9 shows the effect of aspect ratio of the plate on postbuckling response curves of 
piezoelectric plate without damage under different electrical loads. As expected, these results show 
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that the buckling loads are increased by increasing the aspect ratio of the plate. It can also be found 
that the effect of control voltage is more pronounced for the square plate than for the rectangular plate. 

Figure 10 shows the effect of external loads on postbuckling response curves of piezoelectric plate 
with damage and initial deflection. The electrical load is taken as V e = 0 . It can be concluded that the 
larger the external compressive loads, the quicker the development of damage and the more obvious 
the effect of damage on the postbuckling deflection. From Figure 11, the conclusion that the change of 
initial deflection has a small effect on the postbuckled deflection and the damage effect varies slightly 
can be drawn. 

Figure 9. Effect of aspect ratio on postbuckling response curves of piezoelectric plate 
without damage under different electrical loads. 




Figure 10. Effect of external loads on postbuckling response curves of piezoelectric plate 
with damage and initial deflection. 
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Figure 11. Effect of initial deflections on the postbuckling response curves of piezoelectric 
plate with damage. 
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Figure 12 shows the effect of electrical loads on the postbuckling response curves of piezoelectric 
plate with damage and initial deflection. It can be seen that the control voltage has a notable effect on 
the postbuckled deflection and the effect of damage varies greatly. The negative control voltage results 
in a smaller rate of damage development of the plate than that of the same plate without electrical 
loads, so it can be concluded that the positive control voltage can increase the rate of degradation of 
the stiffness of the piezoelectric plate. This can be explained as follows: the effect of positive control 
voltage on the damaged piezoelectric plate is equivalent to a tensile force acting in the in-plane 
direction of the plate, resulting in the acceleration of the piezoelectric plate degradation process. 
Therefore, it demonstrates a prominent effect on the deformation of piezoelectric plate with damage 
under the same loading capacity. 

Figure 12. Effect of electrical loads on postbuckling response curves of piezoelectric plate 
with damage and initial deflection. 
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Figure 13 shows the effect of aspect ratio on postbuckling response curves of the piezoelectric plate 
with damage and initial deflections and Figure 14 shows the effect of thick-span ratio on postbuckling 
response curves of the piezoelectric plate with damage and initial deflections. The electrical load is 
taken as V e = 0 . From Figure 13, it can be seen that the bigger the aspect ratio of the plate, the smaller 
the postbuckled deflection of the plate and the smaller the effect of damage. Figure 14 indicates that 
the larger the thick-span ratio, the smaller the postbuckled deflection of the plate and the smaller the 
effect of damage. 

Figure 13. Effect of aspect ratio on postbuckling response curves of piezoelectric plate 
with damage and initial deflection. 
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Figure 14. Effect of thick-span ratio on postbuckling response curves of piezoelectric plate 
with damage and initial deflection. 
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5. Conclusions 

This paper presents an approach to investigate the postbuckling analysis of piezoelectric plates 
including damage effects using Talreja's tensor valued internal state damage variables and the 
Kachanvo damage evolution equation. The effects of applied voltage, plate aspect ratio, thick-span 
ratio, damage as well as initial geometric deflections on the postbuckling behaviors of the piezoelectric 
plate are investigated. Numerical results show that the nonlinearity of structure has a great influence on 
the postbuckling paths of the piezoelectric plate. The negative control voltage results in the increase of 
the buckling loads and the decrease of postbuckled deflections under the same in-plane compressive 
loads, whereas the positive control voltage decreases the buckling loads and induces larger 
postbuckled deflections. The buckling loads increase with the increase of the thick-span ratio of the 
plate, and the control voltage has a small effect on the postbuckling behaviors of the plate with lower 
thick-span ratio. When the damage and damage evolution are considered, the postbuckled deflection of 
the plate will gradually grow with the increase of the time until the damage reaches a characteristic 
damage state. The external in-plane compressive loads and the applied control voltage have great 
effects on the postbuckled deflections of the plate and the damage development. The negative control 
voltage can decrease the degradation rate of the stiffness of the piezoelectric structures and will 
provide a control mean for the damaged smart structures. 
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